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ABSTRACT 

At the beginning of the reionization epoch, radiation sources produce fluctuations in the redshifted 21 -cm 
background. We show that different types of sources (such as miniquasars, Pop II and III stars, supernovae, 
etc.) produce distinct signatures in the 21 -cm signal radial profiles and statistical fluctuations, through which 
they can be identified. Further, we show that the 21 -cm signal from X-ray emitting sources is much easier to 
observe than was expected, due to a previously neglected pumping mechanism. 
Subject headings: cosmology: theory - early universe - galaxies: formation 



1. INTRODUCTION 



As shown by WMAP observations ( Sper gel et al. 2 006). 
within a half billion years after the Big Bang the Universe 
must have given birth to numerous radiation sources, that by 
z ~ 10 reionized most of the gas in the intergalactic medium 
(IGM). However, the nature of these sources remains a puz- 
zle. Miniquasars, Pop II and III stars, supernovae, and many 
other s have been suggested as the m ain engines of reioniza- 
tion ie.s. Choudhurv & Ferrara 20061) . Unfortunately, direct 
observations may not be able to distinguish between the dif- 
ferent types. For most, the low-energy part of the spectrum 
has the same thermal shape, J{v) oc v 2 , while the upper part, 
hv > 10.2 eV is absorbed or scattered by the neutral IGM. 
However, as we show in this paper, it is still possible to mea- 
sure the shape of the upper part of the spectrum, from the 
source imprint on the 21 -cm signal of the surrounding neutral 
IGM. 

When at redshift z the hydrogen spin temperature, T s , 
is different from the temperature of the cosmic microwave 
background (CMB), 7cmb, the CMB brightness temperature 
presently observed at the wavelength 21(1 +z) cm changes by 
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where S is the local overdensity and x is the ionized hydrogen 
fraction. The spin temperature is determined by (Field 1958) 
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where 7^ is the kinetic temperature, y c is a constant propor- 
tional to the collisional excitation rate (see Zygelman (2005) 
and Liszt (2001) for the contribution of neutral atoms and 
electrons to the collisional pumping) and y a is given by 
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where = 0.416 is the oscillator strength of the Lya tran- 
sition, Ajo is the spontaneous emission coefficient of the 
hyperfine transition, /n is the radiation intensity at the far 
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red wing of the Lya resonance (where the intensity of UV 
photons is unaffected by scatterings) and S a is a correc- 
tion factor, which for th e unperturbed IGM is given by 
dChuzhov & Shapiro 20051) 
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Several authors have previously explored the ways T s can be 
decoupled from 7cmb by the first radiation sources. Nusser 
(2005) has pointed out that, since charged particles have a 
much higher cross-section for collisionally de-exciting the hy- 
perfine transition, T s can be decoupled from 7cmb in moder- 
ately overdense partially ionized regions (which may be pro- 
duced by an X-ray source or following the recombination of 
a fully ionized HII region). In this way a miniquasar with 
luminosity L = 10 41 erg • s" 1 can d ecouple the gas w ithin a 
few hundred comoving kiloparsec ( iKuhlen et al. 20051) . Cen 
(2005) considered a massive starburst in the first galaxies that 
produce a strong UV background. In his model T s is de- 
coupled from 7cmb by the photons slightly on the blue side 
of the Lya resonance that redshift into the resonance within 
some distance (typically up to a few comoving Mpc) from 
the source. However, in this paper we show that in all of the 
above systems the spin temperature is dominated by locally 
emitted Lya photons that can be produced in three different 
ways (collisional excitation, recombination, or cascade). In- 
cluding these photons in the calculation not only boosts the 
intensity of the signal (by orders of magnitude for the X-ray 
sources), but makes a qualitative change in its radial profile. 

In §2, 3 and 4 we recalculate the 21 cm signal from the 
systems discussed above (X-ray and UV sources, and recom- 
bining HII regions). In §5 and §6 we consider the signal pro- 
duced by galactic winds and short-lived radiation sources, re- 
spectively. In §7 we summarize our results. 

2. X-RAY SOURCES 

After traveling an average comoving distance of \(v) rj 
0.1 Mpc ((1 +z)/20)~ 2 (hv/0A keV) 3 through the neutral gas, 
an X-ray photon gets absorbed by an atom. A high energy 
electron produced by the absorption scatters with other atoms, 
electrons, and protons, thereby depositing its energy into heat, 
secondary ionizations, and atomic excitations. Shull & Van 
Steenberg (1985) showed that the division of energy depends 
on the ionized fraction x, with <\>n 
total going to heat, 0.39(1— X 
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On the basis of these results we find that when the initial 
ionized fraction is low (as it is for the primordial gas after 
Z ~ 1000) and temperature is below 10 4 K, the X-rays increase 
7k and x by 
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where e is the average energy deposited per atom. 

The heating rate and intensity of Lya photons produced by 
deexcitations of hydrogen atoms, at a distance R from the X- 
ray source, is given by 
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where L u is the source luminosity at frequency v and is the 
hydrogen atom number density. 

Assuming that the spectrum of high-redshift X-ray sources 
is similar to "ultraluminous" X-ra y sources in nearby gala xies, 
L v w Lois' 1 for 0. 1 < hv < 2 keV (le. g. Miller et al. 20031) . we 
can transform eqs. |6]and0into 
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where f is given in Myrs. 

Figure 1 a shows the evolution of 2 1 -cm brightness tempera- 
ture. Close to the source, where the gas is heated above Tqmb, 
Tt, is positive. Farther out the gas becomes colder, so that Tt 
becomes negative. Still farther out, the pumping radiation be- 
comes too weak to decouple T s from T CMB , and T/, approaches 
a small limiting value that is given by collisional coupling 
only. 

3. UV SOURCES 

The decoupling of the spin temperature may also be done 
by radiation sources that emit a significant fraction of their en- 
ergy in the range between Lya and Ly-limit frequencies. Un- 
like ionizing photons, photons in this range can travel freely 
through neutral hydrogen, until they are redshifted into one 
of the atomic resonances. With the exception of photons that 
are very close to the Ly-limit, which may skip one or more 
resonances, most of these photons are first scattered when 
the closest resonance is reached. Thus, only photons emit- 
ted between Lya and Ly/3 frequencies reach the Lya reso- 
nance, which they leave only after ~ 10 6 scatterings. By 
contrast, photons emitted above the Ly/3 frequency typically 
scatter only a few times before being destroyed by cascading 
to lower levels. A fraction of resonant photons destroyed by 
cascade (0 for Ly/3, 0.26 for Ly£ and ~ 1/3 for higher res- 
onanc es) produce Lya photons as one of the c ascade prod- 
ucts (iHirata 20061 iPritchard & Furlanetto 20061) . Since the 
frequency range between the Ly7 and Ly-limit is 2.2 times 
smaller than that between Lya and Ly/3, Lya photons pro- 
duced directly through redshifting are generally more numer- 
ous than those produced by cascade. However, the latter pho- 
tons are typically produced within a much shorter distance 
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FIG. 1. — The brightness temperature versus comoving radius around a 
radiation source at z = 20. a) X-ray source with L = 10 42 erg/ s and the age of 
t = 1 (solid), 2 (dash-dotted) and 4 Myr (dotted line), b) UV source with total 
luminosity of, respectively, L = 5 X 10 41 (dashed), 5 X 10 42 (dash-dotted) and 
5 X 10 43 erg/s (solid line) between Lya and Ly-limit frequencies, c) Galactic 
wind with L a = 10 41 (solid), 10 40 (dashed) and 10 39 erg/s (dash-dotted line). 



from the source (on average the photon between n and n + 1 
resonance travels a distance of ~ «" 3 //"'c/2 before reach- 
ing a resonance), which makes them a dominant pumping 
mechanism within a few Mpc from the source. Including 
these photons in the calculation not only boosts the pump- 
ing radiation intensity, but also changes its radial profile. 
Whereas normally the radiation intensity is proportional to 
R~ 2 , the cascades change this dependence to approximately 
R~ 7 ^. In addition, they also result in the discontinuous jumps 
in pumping radiation intensity and 21 cm brightness signal at 
R n = (n~ 2 -(n+ l)~ 2 )H~ l c. Figure lb shows the radial profile 
of 21 -cm brightness temperature of the unheated IGM around 
UV sources. In principle, the resonant scatterings affect the 
gas temperature, which therefore becomes time- and position- 
dependent (Chuzhov & Shapi ro 20051) . However, unlike the 
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FIG. 2. — The brightness temperature versus comoving radius of the re- 
combining HII region that was ionized at z = 20. The solid, dashed and dash- 
dotted lines correspond to z = 19, 17 and 15. The numbers above and below 
the lines show, respectively, 7j./ 1000 K and x at a given redshift. 
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X-ray photons, the UV resonance photons generally gain or 
lose only a very small fraction of their energy in exchange 
with the gas, so in Figure lb we neglect this energy transfer. 

4. RECOMBINING HII REGIONS 

The low-mass galaxies that are first to form in the standard 
ACDM model are expected to be very vulnerable to feed- 
back. During brief and intensive periods of star formation, 
such galaxies may be completely disrupted by SN explosions 
or photoevaporation. Furthermore, reionization may also de- 
lay gas accretion and star formation in the neigbouring halos. 
Consequently, the H II regions produced by the first galaxies 
may be given a pause, allowing them to cool and recombine. 

Unless the gas is extremely overdense, its thermal evolution 
as it recombines is determined mainly by inverse Compton 
cooling and radiative recombination. Right after ionization 
the temperature of the gas in the IGM depends on the shape 
of the ionizing spectrum. If the shape changes appreciably 
during the life-time of the ionizing source or if the source life- 
time is smaller or comparable to the gas cooling or recombi- 
nation times, this would be reflected in the radial profile of the 
gas temperature and ionized fraction. Also, a strong jump in 
both temperature and ionized fraction can be produced across 
a He II ionization front. Even after the radiation source is 
dead, the difference in the ionized fraction inside and outside 
a former He III region would persist for quite a long time. 

Figure |3 shows the evolution of the 21-cm signal from 
the recombining H II region, produced at z = 20 by a short- 
lived source whose spectrum is similar to a black body with 
a surface temperature of 90,000 K. Initially the gas bright- 
ness grows almost linearly with the hydrogen neutral fraction. 
Subsequently, when the ionized fraction falls below ~ 0.2, 
the intensity of Lya photons (which dominate over collisional 
coupling for x above ~ 0.05) produced by recombinations is 
insufficient to fully decouple T s from Tqmb, and, unless addi- 
tional Lya photons are provided by exterior radiation sources, 
the signal may drop. 

5. GALACTIC WINDS 

Even though large numbers of Lya photons may be pro- 
duced by the overdense gas in galaxies via recombination or 
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FIG. 3. — The brightness temperature versus comoving radius around a 
short-lived UV source at z = 20. The intensity of the source is assumed to 
vary as L = 5 X 10 43 erg/s X sin(7rt/3M>r) for t < 3 Myr. The solid, dashed 
and dash-dotted lines correspond to t = 3.5, 4 and 5 Myr. 




FIG. 4. — Power spectrum of the 21 cm signal at z = 20. The pumping 
radiation is assumed to be produced by X-ray sources in halos with T yII > 10 4 
K or r v j, > 5 ■ 10 3 K (solid and dotted lines), or by Pop III stars in halos with 
Tyir > 10 4 K or r v j, > 5 ■ 10 3 K (dashed and dash-dotted lines). We assumed 
>' Q < 5. 



excitation of hydrogen atoms, these photons do not typically 
affect the 21 cm emission in the neighbouring IGM, as they 
redshift out of resonance close to their origin. This limita- 
tion, however, may be bypassed if the peculiar velocity of the 
gas is very large. Lya photons emitted by a medium mov- 
ing with peculiar velocity V may reenter the resonance within 
R = V /H from the point of origin. Large peculiar velocities 
involving a large mass of gas may indee d be generated by one 
or more supernovae explosions (le.g . Ostriker & Cowie 198 It 
iKitavama & Yo shida 2005; Wilm an et al. 20051) . 

A wind expanding with velocity V s that exceeds 100 km/s 
shockheats and compresses the surrounding gas into a thin 
shell. The cooling in the outer parts of the shell, where the 
temperature is above ~ 10 6 K, is dominated by inverse Comp- 
ton scatterings and bremsstrahlung. In the inner parts, where 
the temperature drops below ~ 10 5 K, most of the energy is 
radiated away in Lya photons. By considering a thin shell 
expanding with velocity V s at a radius R s , we find the distribu- 
tion of the photons that are at resonance in the local gas frame 
to be 

Figure lc shows the 21 cm signal produced by a galactic wind 
for different radiation intensities. 

Apart from supernovae, thin shells expanding with large pe- 
culiar velocities can be produced by large voids forming from 
steep initial density perturbations. The distribution of Lya 
photons produced by such shells can still be described by 
equation ITOl However, in this case the ratio V S /HR S (which 
for shells driven by s upernovae can be very larg e) is typically 
between 1.2 and 1.5 (Chuzhov & Nusser 2002). 

6. SHORT-LIVED RADIATION SOURCES 

In the examples given in the previous sections, we have 
neglected the possible time-dependence of the luminosity of 
radiation sources. However, if the luminosity changes over 
a relatively short time-scale, this will be reflected in the 21- 
cm signal. Since the speed of light is finite, regions that are 
observed simultaneously are pumped by photons emitted at 
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different times. Therefore, by measuring the 21-cm signal 
from different regions around the source, it is possible to re- 
construct the evolution of the source spectrum and intensity. 
Figure |5] illustrates the evolution of the 21-cm profile after a 
starburst that lasts a few million years. 

7. SUMMARY 

We have shown that different radiation sources produce dis- 
tinctly different 21-cm profiles, which can be sensitive not 
only to the source spectrum, but also to its history (§2 and 
6). Furthermore, differences can also be identified in the 
combined signal from multiple radiation sources of different 
types. Figure|4]shows the power spectrum of the 21-cm signal 
produced by either the X-ray or the UV sources contributed 
by a biased distribution of halos (the power spectrum calcula- 
tion is based on the formalism developed by Barkana & Loeb 
(2005a)). Predictably, the X-ray sources, which typically in- 
fluence smaller regions, produce greater fluctuations on small 
scales than do Pop III UV sources. This difference between 
X-ray and UV can be even larger if Pop III is replaced with 
Pop II stars, which produce lower proportion of higher n res- 
onance photons (Barkana & Loeb 2005b; Pritchard & Furlan- 
etto 2006). 

Assuming that at z ~ 20 the typical distance between the 



first radiation sources is a few Mpc, their individual detec- 
tion requires an angular resolution of about one arcminute (or 
baseline of ~ 15 km) and frequency resolution of ~ 100 kHz. 
For a beamwidth and bandwidth just small enough to resolve 
such objects, a detection of an X-ray source with L ~ 10 42 
erg/s (which equals the luminosity of one or more black holes 
of total mass ~ 10 4 M Q shining with Eddington luminosity) 
requires a sensitivity of A/T ~ 10 3 m 2 K -1 / y/ '(f bs/ 'I0 3 hours), 
where f b s is the integration time. Thus it appears that SKA, 
whose planned sensitivity is a few times 10 3 m 2 K _1 at the rel- 
evant frequencies, would be able to detect such X-ray sources 
individually. The detection of the fluctuating 21-cm back- 
ground produced by the superposition of the first sources is 
even more promising. For example, the curves in Figure |4] 
can be easily distinguished at scales k ~ 10~ 2 Mpc" 1 by LO- 
FAR, SKA, or PAST for 4 weeks of integration time (see Zal- 
darriaga, Furlanetto, & Hernquist 2004 for details on power 
spectrum error estimation). 
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